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We study the current, the curvature of levels, and the finite temperature charge stifTness, D{T, L), 
in the strongly correlated limit, U ^ t, for Hubbard rings of L sites, with U the on-site Coulomb 
repulsion and t the hopping integral. Our study is done for finite-size systems and any band filling. 
Up to order t we derive our results following two independent approaches, namely, using the solution 
provided by the Bethe ansatz and the solution provided by an algebraic method, where the electronic 
operators are represented in a slave-fermion picture. We find that, in the U — oo case, the finite- 
temperature charge stiffness is finite for electronic densities, n, smaller than one. These results 
are essencially those of spinless fermions in a lattice of size L, apart from small corrections coming 
from a statistical fiux, due to the spin degrees of freedom. Up to order t, the Mott-Ifubbard gap 
is Amh = U — 4t, and we find that D{T) is finite for n < 1, but is zero at half-filling. This result 
comes from the effective flux felt by the holon excitations, which, due to the presence of doubly 
occupied sites, is renormalized to ^'^■^■^ — (f>{Nh — Nii)/{Nd + N^), and which is zero at half-filling, 
with Nd and Nh being the number of doubly occupied and empty lattice sites, respectively. Further, 
for half-filling, the current transported by any eigenstate of the system is zero and, therefore, D{T) 
is also zero. 
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I. INTRODUCTION 



The transport properties of strongly correlated elec- 
trons in low-dimensional conductors has been a subject 
of experimental and theoretical interest for over twenty 
years. Low-dimensional conductors show large devia- 
tions in their transport properties from the usual single- 
particle description. This suggests that electrooic-corre- 
lations play an important role in these systemsETLl, even 
if the correlations are smallQ. Solvable one-dimensional 
many-electron models such as the Hubbard chain are of- 
ten used as an approximation for the study of the proper- 
ties of quasi-one-dimensional conductors. Although the 
Hubbard chain has been diagonalized long agdj'B, the in- 
volved form of the Bethe-ansatz (BA) wave function has 
prevented the calculation of dynamic response functions, 
these including the charge-charge and spin-spin response 
functions and their associated conductivity spectra. In- 
formation on low-energy expressions for correlation func- 
tions can be-obtained by combining BA with conformal- 
field theorvEl. On the ojtlier hand, several— approaches 
using perturbation theoryllil, bosonizationll3ll3, the pseu- 



doparticle fcpmalism&lli, scaling methodsS, and spin- 
wave theorylia have been used to investigate the low- 
energy transport properties of the model away, from half 
filling and at the metal - insulator transitionB. Limited 
information on the transport properties ajtiinite energies 
has been obtained by numerical methodst£l 

Recently, a series of comparative numerical and ana- 
lytical studies have explored the differences in the trans- 
port priipeiiies between integrable and non-integrable ID 
modelscJcj, at finite temperature. Most of these stud- 
ies have dealt with generalizations to finite-temperati|tce 
of Kohn's zero-temperature concepts and approachEj. 
Zotos and Prelovsekcj have introduced the concepts of 
ideal insulator and ideal conductor at finite teurpeTeLtures. 
These concepts refer to the temperature dependence of 
the real part of the optical conductivity ar{^,T), which 
is given by 



ar{uj,T) = 2ttD{T)5{uj) 



reg 



(1) 



where we have taken ?i = = 1, e is the electron 
charge, cu is the frequency of the electric field, and T 
is the temperature. The quantity D{T) is the charge 



1 



stiffness and characterizes the response of the system to 
a static electric field, within linear response theory. Ac- 
cording to Kohn's zero-temperature criterion, the value 
of D{Q) can be used to distinguish between an ideal in- 
sulator -D{0) = - and an ideal conductor - D{0) ^ 0. 
The quantity areg{'-^,T) is the regular part of the con- 
ductivity and describes the absorption of light of finite 
frequency uj by the system. 

In this work we are concerned with the one-dimensional 
Hubbard model. The response of the energy eigenvalues, 
Em, of the model to an external flux (/) = ■AL {A is the 
vector potencial along the chain) piercing the ring can be 
used to determine D{T). 

The charge stiffness D{T,L), for a finite-size sys- 
tem of l^AStt L, can be evaluated as a thermodynamic 



quantitj 



and 



and is given by 



(fE„ 



d{<p/Lf 



(2) 



(3) 



i=0 



where Z is the partition function and Dm is the curvature 
of the eigenenergy £',„. This equation is derived using 
finite temperature linear response, and was obtained for 
the first time by Castella, Zotos, and Prelovsek in connec- 
tion with the problem of diffuspe or ballistic transport 
of a particle in a fermionic bathc3. 

Another equation for the charge stiffness D{T) can be 
derived starting with the partition function Z{(j)). Taking 
the second derivative of the free energy !F = — TlnZ(0) 
in order to 0/L we arrive to an alternative relation for 
D(T) 



D{T,L)^Di{T,L) + D2{T,L) 



(4) 



where 



1 d^T 



0=0 



and 



Jr. 



d{<j>/L) 



(5) 



b=0 



The terms Di and D2 represent the thermodynamic and 
the current contributions, respectively. The thermody- 
namic contphution, Di{T), vanishes in the thermody- 
namic limitcj. |— I 

Zotos and Prelovsekcil, based on numerical work done 
in small sistems, conjectured, for intcgrable systems - 
namely the spinless-fermion and the Hubbard models -, 
that D{T) should be zero if D{0) is also zero, but finite 



otherwise. Very recently, ZotoscSI has shown that this is 
so for the isotropic Heisenberg model, and some of the 
present authors have showncj that Zotos conjecture also 
holds for the t—V model, in the strong interacting limit 
V ^ t, with V the nearest-neighbour Coulomb integral. 

Very recently, Kirchncr, Evertz, and Hanke have stud.- 
ied D{T, L) using Quantum Monte Carlo simulationsEil. 
They have found that pome conjectures first proposed 
by Zotos and PrelovsekEj, regarding the issue of ideal- 
insulator behavior in one-dimensional integrable systems, 
should not be true, in quali-tative agreement with a ther- 
modynamic Bethe-ansjatzL^ calculation performed by Fu- 
jimoto and KawakamiE3. 

At finite values of L and at zero temperature, the full 
calculation of -D(0, L), for any value-of U ^ was first consid- 
ered by Stafford, Millis and ShastryllZl. These authors de- 
rived scaling relations for 15(0, L). In particular, at half- 
filling and for the strong coupling regime, [/ — > cx), we 
have £1(0, L) ~ {-l)^/'^+^L^/'^AUeyi]i[-L/£,[U)], where 
A is a constant. The study of persistent currents jo{4>/L) 
in the ground state for finite size Hubbard rings and its 
relation with mesoscopic transport in arrays of quantum 
dots has been cairsidered by KusmartseVj^eisz, Kishore, 
and TakahashiEJ, and by Yu and FowlertS. 

It is interesting to remark that, as in a Fermi liquid, 
in the thermodynamic limit, D{0, 00) can also be de- 
rived by nieans of kinetic equations for the elementary 
excitationslla. 

In this paper we present an analytical study of jm. 
Dm, and D{T,L) for Hubbard rings of different sizes. 
Our results refer to infinite and very large values of U. 
In addition to the BA based cepresentation, we also use 
a slave-fermion representationEZI, both of them leading to 
the same physical results. 

The paper is organized as follows. In Sec. |l|we present 
the solution of the BA equations with twisted bound- 
ary conditions. We use these equations to compute the 
mean value of the current operator in any eigenstate of 
the model and the corresponding curvature of the en- 
ergy level. These results allow us to determine the finite- 
temperature behavior of the charge stiffness. We do our 
calcutions in the U = 00 and U ^ t cases. In Sec. 
HI we rederive the results of Sec 0, using an algebraic 
approach, where the electronic operators are described 
in a slave-fermion representation, and the same physi- 



IV 



cal quantities are computed up to order t. In Sec 
we make a comparison between two charge gapped sys- 
tems, at half-filling: the Mott-Hubbard insulator (repre- 
sented by the Hubbard model) and the band insulator 
(represented by a dimerized lattice model). In Sec. ^ we 
present our conclusions. 



II. CHARGE STIFFNESS D{T) FOR A 
FINITE-SIZE SYSTEM 
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A. Bethe-ansatz equations 

In this section we study the charge stiffness at finite 
temperature for one-dimensional Hubbard rings of finite- 
size L. The Hamiltonian for the Hubbard model is given 

by 



H 



H.c] 



(6) 



where c|^(cjcr) creates (annihilates) one electron with 
spin a (here and when used as operator index, a =t: i? 
and a — ±1 otherwise), fij^a — Cja'^ja- is the number 
operator at site j, and cl+io- = cio-- 

In the presence of a time-dependent vector potencial 
A{t) the hopping amplitude t changes according to the 
well-known Peierls substitution 

U+i,,{A{t)) = ti+i^, exp[-ie(i?,+i - Ri)A/nc] , (7) 

where Ri is the position of the site i on the lattice. It has 
been possible to solve the Hamiltonian (^) with the ad- 
ditional hopping phase exp([— ie(i?i+i — Ri)A/hc\). For 
convenience we write A = Tictj) / {aeL) , where the lattice 
spacing a is given by a = -R^+i — Ri- More generaly, 
a spin-dependent Aa vector potencial can be introduced 
and the model (|^) can be solved by means of the co- 
ordinate BA either with twisted or toroidal boundary 
conditior*. both approaches giving essentially the same 
resulta23'E2l. One obtains the energy spectrum of the 
model parameterized by a set of numbers {kj, A^} which 
are solution of the BA interaction equations given by 



N 



and 



TT sinfcj 
sinfcj 



Aa + iU/4. 



As - iU/A 



,N). 



(8) 



N 



TT sinfcj - As + iU/4: 
JJ^ sin kj - As - ill/ 4 



N 



ji Ap-As + tU/2 ^ ^ 

J-J-^ A^ - A,- - tU/2 ' ^ 



(9) 



The above equations have both real and complex so- 
lutions for the rapidities kj and.. Ag. However, previous 
studies of the (j)^ ^ problen£3o have only considered 
the real BA rapidities solutions of Eqs. (^) and (||) which 
refer to the low energy excitation spectra. The general 
solution of Eqs. (||) and (^ for a finite system. a£ size L 
has been obtained previously in the literaturerHP'^l and is 
given by 



kjL ^ 2t:I!I + - X! X! ^ ^^'^ 
7 i'=i 

-EE^tan 

7>0j' = l 



_i f sinkj/u — Rs.-y,j 



7 + 1 



^ ( smkj/u — Rc.-i,j' 



(10) 



2L sin 



(^^1^ + (Ren.] + 1/")^ 



E^" ^ 1 / sin kji lu — Rc-y 1 
2tan"^ ' ^— 



and 



^ 2 tan" 



'E - Rc.,-i',r) 1 (11) 

7'>o i'=i 



_l ( Rs,-y.j — smkjt /u 



i'=i 



1 + 7 



2^/;^^ + E E ^-,+in'+iiRsn,3 - Rsn',r) ■ (12) 

7' J'=l 

In the above equations we have considered the case 
4>a = <f)^a and u — U/{4t). 

The functions Q-y^-yr{x) [and Q-y^i^^i{x)] of Eqs. 
( pO| ) , (O) , and (|l^) are defined in Ref. [l^, and are sums 
of tan (a;). The following definitions for the real part 
of the rapidities, AJJ+^/u = Rs.-yj (with n + 1 = 7 and 
a — j), AJ^/u — Rc.-y,j (with n = 7 and a = j), and 
7 = 1, 2, . . . , CX3 for the sums and 7 = 0, 1, 2, . . . , cx) 
for the Ns^-y sums, allow us to recover Takahashi's formu- 
lae for (j) — (n. Here and often below we use the notation 
c = c, 0, which allows the c, 7 sums to run over 1,2,3,... 

The relevant numbers 1^, Ij''^: and /|''^ which appear 
in going from Eqs. (|) and (|) to Eqs. (0), (|ll|), and (|l|) 
are the quantum numbers which describe the Hamilto- 
nian eigenstates. In Table | we give a classification of the 
type of excitations described by these quantum numberi 

The numbers I"''' can be integers or half-odd integersi 
if the numbers are even or odd, respectively, where 



7—0 7—1 

Ncn = 1 + L - N + Nc-y , 



(13) 



The spacing between adjacent quantum numbers I^''^ , 
is always one, and independent of the value of the in- 
teraction U. It is therefor e n atural to interpret qj''^ — 
2ttI"''^ /L as a momentumcS, and the rapidities Ra,j,j 
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and kj as functions of q"''^ ■ The total number of a elec- 
trons, N^, is given by the constraints 

iVx = A, + A, , (14) 

and 

N = Nc + 2A,, (15) 

where 



7=1 

^ = ^(7 + l)iV.,7, 

7=0 



(16) 



and the numbers Ij, Ij''^, and Ij''^ belong to the intervals 

|/;i<|, 
\i;'''\<^{L-N + 2A,~Tj), 

where TJ (with a = c, s) are given by 

7' = 1 
7'=0 



(17) 



(18) 



with f^ .^, = 2min(7,7') - S^^-y> and i^ .^, = 2min(7 + 
1,7' + i)- (5^,7'. 

All the eigenstates considered in the BA solution of the 
model are described by the different occupancies of the 
quantum immbers If''' ■ For example, the ground state is 
describcdtS by a compact symmetric occupancy around 
the origin of the numbers Ij and Ij'^, and by zero oc- 
cupancy for the numbers I^"' and The general 
situation is given in Table 

Tiie energy and momentum eigenvalues are given 
byi0 



E{L,(j),U,N„) = -2t^coskj 

Na,-, 



7=1 j = l 



(19) 



and 



L 3 ^ ^ L 

j = l 7 j = l 7 j=l 

^(A + 7r2^7Vc,7- 

7 



(20) 



It is interesting to remark that the flux contribution 
for the total momentum P differs from its contribution in 
the independent electronic problem, for which we would 
obtain a contribution of the form P^ = (t)N/{2L). The 
general solution of Eqs. (10)-(^2|) for arbitrary system 
size L, electron numbers Na, and Coulomb interaction 
[/ is a very difficult problem. Close to and at half-filling 
Stafford et al. computed the-^ero temperature charge 
stiffness for a system of size ZJlZl. To compute the charge 
stiffness at finite temperature we need to compute the 
mean value of the current operator or the curvature of 
levels of all the eigenstates of the model and not only 
that of the ground state. Below, we solve the problem 
explicitly for a system of size L in the limits U = oo and 
U ^ t. For simplicity we shall consider the case where 
the quantum numbers /"'''' are symmetrically distributed 
arround the origin. 



B. The U = oo case 

At J7 = cxD the charge strings have infinite energy and, 
therefore, drop out of the problem. In physical terms this 
means that the states with one or more doubly occupied 
sites are not permited. The rapidites i?s,7,j decouple 
from the charge degrees of freedom and all the spin ex- 
citations are degenerateEJO'Ej. The s,0 spinon and the 
s, 7 > spin-string excitations have a flat dispersion re- 
lation. The equations (p^)-(|l^) reduce to 

k°°L = 27r/| + + ^ £ 2 tan-i (R^.r/il + 1)) 

7 j' = l 

+ 5]g2tan-i(i?-_^.,/7) , (21) 

7 j' = l 



(L-iV,)2tan-i(i?-^^./7) = 
2^If^ _ 270 + ^ ^ - Ki',r ) ' (22) 

7' j' = l 



and 



iV,2tan-i(i?- ^^7(1 + 7)) 



= 2^j;'^ + E E ©7+1.7'+! (^^7,. - ) ■ (23) 

7' = l 

We see from the structure of Eqs. (pll)-(p3|) that the 
Hubbard model Hilbert space decouples into a product 
of three Hilbert subspaces, each one associated with a 
different Hamiltonian. The latter are a chain of length 
L and Nc spinless fermions, an Heisenberg spin one-half 
chain of length Nc , and an Heisenberg spin one-half chain 
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of length L — Nc- This is, however, a very dehcated de- 
couphng in what regards the thermodynamic properties. 
Since in a thermodynamic calculation N^. must vary be- 
tween Nc = L and Nc — 0, the length of the two Heisen- 
berg chains also varies. 

The energy eigenvalues are given by 



E{L,(j),N„) = = -2t y cosfc 



(24) 



where the subscript m labels a given eigenstate of the 
model. Similar equations to Eqs. (pTl)-(p4[) have been 
derived for the ground statei in the study of persistent 
currents in finite-size ringscil. If we use Eqs. ( p2| ) and 
( p3| ) in Eq. (|2^) we obtain a solution for fc°° in terms 
of the quantum numbers ^ and only. This solution 
reads 



,00 ^ 2vr £-iVc-E,27iVc., 
^ L ^ L L~Nc 



271 

LMc 



7 j=i 



2n 



L{L - Nc) 



(25) 



fermions. Small differences are to be expected due to the 
statistical flux term tAf*"* 



2tt 

LNr 



EE 

7 i=i 



ir. 



(29) 



If we consider the special case of a fully polarized sys- 
tem the statistical flux term vanishes and the system is 
equivalent to free spinless fermions in a ring of size L. 

To illustrate the differences between free spinless 
fermions and the Hubbard model at J7 = 00 when 
^stat we plot in Figure D{T,L) for the C/ = c» 
Hubbard and for the spinless free-fermion models in a 
ring of L = 12 for N^ = Ni = 3. For this case we have, in 
the Hubbard model, two types of spin-string excitations, 
those with 7 = 1 and 7 = 2. The details of the states 



available are given in Table III. The results of D{T, L) in 
Figure ^ are normalized by the zero-temperature charge 
stiffness 15(0, L) which, for TV even, is given for the Hub- 
bard model by 



7 j=i 



sin(7r7V/L) 1 
sm(n/L) L 

and for free spinless fermions by 



(30) 



Equation (|2^) shows that the spin degrees of freedom are 
still coupled to the charge excitations through the quan- 
tum numbers ■ These act as a ficticious flux pierc- 
ing the ring of spinless fermions. The fact the charge 
and spin degrees of freedom are not completely decou- 
pled introduces interesting statistical consequences in the 
classificatiop-joLthe excitations according to the Haldane 
criteriumyee 

The eigenstates with non-zero occupancy of the num- 
bers Nc.-f have infinite energy (for U — 00) and therefore 
they drop out of the problem. That is, k°° has the sim- 
pler form 



L J 



(j) 27r 
L ^ Tn, 



■EE 

7 j=i 



(26) 



The curvature of levels D„i{L) for an eigenstate m is 
given by 



1 d^E„ 



2id(0/L)2 L^j^^ 



7E'^(^ 



(27) 



where the sum is over the configuration of quantum num- 
bers Ij defining the eigenstate m and d{Ij) is given by 



d{n) = t cos kj 



t cos k° 



d ki 



+ t sin /c, 



d'^kj 

WW) 



(28) 



From Eq. (E6f) we expect the behavior of the charge stiff- 
ness D{T, Ufto be essencialy the same as for free spinless 



Do[Q,L) = Du{Q,L) cos-, (C/ = 0) . (31) 

From Eqs. (^^ and (|2^) it is clear that at half-filling the 
charge stiffness D{T) is zero at any temperature, since 
the only available state is the ground state (apart from 
a large massive spin degeneracy). In this limit the spin 
degrees of freedom s, 7 introduce a large degeneracy but 
they do not contribute to the transport of charge. 



FIG. 1. The charge stiffness D{T, 12)/D(0, 12) for L = 12, 
7V| ~ —3. The result for free spinless fermions is the 
dotted line, the full line is the result for the U = 00 Hubbard 
model. 
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We see from Eq. (Q) that the charge stiffness can be 
obtained from a thermodynamic contribution plus a con- 
tribution from the thermal average of the square of the 
mean value of the current operator, jm, that we have 
denoted by D2{T,L). 

In the limit U = oo the mean value of the current 
operator is easily computed in a finite-size system and 
reads 



electrons, whereas their holes, the holons, transport mi- 
nus the charge of an electron. 



E 



2isin/c°°. 



(32) 



In Figure ^ we present the results for D{T,L) and 
D2(T,L), for L = 18 and ^ Ni = 3. We see that 
for T > t/10 D{T,L) and D2{T,L) coincide (T* ~ t/10 
being a crossover temperature for L = 18). On the other 
hand, at low T D2{T,L) is much smaller than D{T,L), 
for finite-size sytems. We have checked that D{T, L) and 
D2{T, L) coincide at lower and lower temperatures, as L 
increases. In the thermodynamic limit the contribution 
Di{T) from the free energy vanishes and D(T) = D2(T). 
Similar results have been obtained by Naroaimy, Millis, 
and Andrei for the spinless fermion modelE3, by exact 
diagonahzation of small clusters (up to L = 14). 

The behavior discussed above can be easily under- 
stood. For a given L, as T increases the contribution 
of the groundstate to the thermodynamic average de- 
creases. Since jo (the current in the groundstate) is zero, 
D2 (T) eventually has to start decreasing to zero at some 
crossover temperature, where the thermodynamic weight 
of the groundstate dominates the Boltzmann average. 
When L increases the number of low energy states above 
the GS also increases. This reduces the thermodynamic 
importance of the jo to the thermodynamic average and 
the crossover temperature T* is shifted to lower values. 

We now show that the results for D2{T, L), away from 
half-filling, for finite-size systems, and low temperature, 
are strongly dependent on the fact that all the spin exci- 
tations have been taken into account. This contrasts with 
the behavior of D{T, L) based on the curvature of levels 
given by Eq. (p8|). In Figure || we plot the contributions 
to D{T) due to the several spin excitations - spinous and 
spin strings - using the curvature of levels and the mean 
value of the square of the current operator. From Figure^ 
we see that if we consider only the spinon s, excitations 
the contribution from these states to D{T^ L) (triangles) 
almost coincides with the full D{T, L) in the whole range 
of temperatures. On the other hand, the contribution 
to D2{T, L) from the same set of states is very small at 
low temperatures. If the s, 1 spin-string excitations are 
included, we see that D2{T, L) is now above D{T, L) for 
low T. As the system size increases the contribution from 
the spin excitations decreases. 

It is interesting to remark that the electric charge Cc 
transported by one c excitation is given by Cc = — 1, at 
U — 00. These c excitations are the c pseudoparticles 
of Refs. p^ , p^ , the holons being their holesEj. Therefore, 
they transport the same electric charge as the original 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 

T/t 

FIG. 2. The charge stiffness, D{T, L)/ D{0, L) (line), and 
the thermal average of (jm)^, D2{T, L)/ D{0, L) (dots), for 
the Hubbard model. The system size is L = 18, the num- 
ber of electrons is Nt; — — 3, and the on-site Coulomb 
repulsion is U — 00. 



FIG. 3. Contributions of the several types of spin exci- 
tations to D{T) and D2{T). The notation is the following: 
D{T, 18) (fuU line); D2{T, 18) (dotted line); D{T, 18) includ- 
ing the s,0 excitations only (triangles); -D2(T, 18) including 
the s, excitations only (squares); D2{T, 18) including the s, 
and s, 1 excitations only (circles). The system size is L = 18, 
the number of electrons is A''| = A''| = 3 and the on-site 
Coulomb repulsion is U — 00. 
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C. The [/ > f case: 0(t) and 0[t^ /U) corrections 



In the subsection fl B| we have considered the case of 
U ^ oo where charge-string excitations are not allowed. 
In this section we consider the case where U ^ t and 
we solve Eqs. (|l|), (|l|), and ® to order t^/U. 

In this case, c, 7-excitations can take place, since their 
energies are no longer infinity. In Table IV we list the 



type of excitations that can take place in the system for 
L = 18 and iV 6. 

In this limit, the Hilbert space of the Hubbard model 
can be decomposed into several subspaces each of energy 
of order jU (7 = 1,2,...) relatively to the ground state. 
In the original lattice this means that we allow for one, 
two, doubly occupied sites. These doubly occupied 
sites are related to the charge-string excitations c, 7. For 
example, in the subspace of energy U we can have only 
one c, 1 string excitation and in the subspace of energy 
2U we can have two c, 1 or one c, 2 string excitations, 
and so on. On the other hand, the conservation of the 
electron numbers in these excitations, imposed by Eqs. 
( |l^ ) and (|l^) , show that the number of occupied /| is re- 
duced by two and by four, respectively. From Eq. (^9|) is 
simple to show that the energy of these string excitations 
is given, to order t, by 

Em = -2t^ cos/cf + UJ2"f^c^'y + 0{f /U) . (33) 
j=i 7 

Equation ( ^3|) has first been derived by WoynarovichE§, 
starling from the ground-state equations of Lieb and 
WuS. In the limit > t and to order /U eqs. (0)-(|l|) 
reduce to 

kjL = 27r/| + + X! X! ^ (-^^,7 J' h) 
7 j'=i 



8t 



-sinfc,^^ 



7 



7 J 



-^^2tan-i(i?,,^,,v/(7 + l)) 
7 j'=i 



8t^ 

'u'' 



7+1 



fr{ (l + 7)2 + (i?.,,,,)^ 



(34) 



(L-iVe)2tan-i(i?,^^,,/7) = 

= 27r/;'^ -2-f(b + J2J2 Q7,7'(^c,7j - Rc,Y,j') 
7' 

and 



iV,2tan-i(i?,,^j/(l-|-7)) 



27r/|'''+^ X! ©7+l:7' + l(-Rs,7,i ~ ^s,7'j') 
7' j' = l 

AT, 



8t 



7 + 1 



C/(l + 7)2 + (i?.,,,,)^^ 



^sinfcj. (36) 



Equations (|3^) , (|35|) , and (^ can be solved explicitly by 
an iteration prpcedure introduced first for the ground- 
state equationstil. The solution, to order t^/U, it is sim- 
ple to obtain and reads 



k, = kr + s^ + s^ 



(37) 



with Sj and Sj given by 



LUN, 



^sinfcj 



LU 



LU{L - N,) ^ 



^sinfcf, (38) 



and 



7 j= 



N, 

EE 

7 j 



1 7^ + (^c°r-y.,)' 



(39) 



(40) 



and given by Eq. (g5|). The sums in Eqs. (^, (|3S|), 
and ( |4(]| ) depend on the distributions of the numbers 1°'''^ , 
that is, on the considered eigenstate. The Bf^ term rep- 
resents the energy of an isotropic Heisenberg chain, of 
length Nc, with the spin-string excitations included. The 
Bm term has the same form as B^ and can be thought 
as an effective Heisenberg chain, of length {L — Nc). At 
odds with Bm, the excitations represented by i?^ have 
a charge gap relatively to the ground state of the sys- 
tem. The smallest gap, A{L,U), is associated with the 
creation of a single c, 1 excitation above the ground state 
and, for a finite-size system, is given by 

A(L, [/) = [/ + 4< coa{Ncn/L - n/L) 
8t2(iVc - 2) 



LU 



-(Bt + l) 



8t\^, , ^^ sin(27rjV,/£-47r/L) 
+ LU^^' + sin(27r/L) 



8t^ 



-B'n N, 



sm{2'KNc/L) 
sin(27r/i) 



oj-2 



(41) 
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where Bq and Bf represent B^^ for the ground state 
and for the excited state considered, respectively. When 
L oo and Nc — > L, we obtain the Mott-Hubbard gap 
Amh = U -U + 8e\n2/U + 0{l/U^). 



The numbers R° 



tion of Eqs. (|22| ) and 
that contains the flux 



and RT,'y,j ^i'^ given by the solu- 
sl), respectively. The term in 
has the form 



L-Nr ~ L ' 



(42) 



Due to the conservation law (^_5|) it is clear that at half 
flUing we have ^""ff /L = and % in Eq. @ only 
depends on 4> through the term. In the asymptotic 
limit U ^ t the sum 7^c,7 is identified with the total 
number of doubly occupied sites Nfj, that is 



(E 



C:7 



+ 0(l/C/2), (43) 



second case, we see that D{0, L) has not corrections of 
the order t'^/U. In fact, the asympotic result shows that 
D{L,U) cx (4.38/L/)^-i. 

We stress that the energy ( ^ only corresponds to the 
true ground state of the model for small values of (j). 
When (j) increases level crossing takes place and there- 
fore the energy (^) does not correspond to the absolute 
ground state of the system. Furthermore, for = 
and L — si— .oo, we recover the results of Carmelo and 
Baeriswylo for the ground-state energy. 

If the corrections ( ^5| ) to the energy are neglected, the 
energy eigenvalues depend on only. We then see, from 
Eq. (Q), that the current dEm/d{(j)/L) and the curva- 
ture of levels d'^ Ejn/d {(j)/ L)'^ are all proportional to z, 
that is 



d Eri 



d(0/L) Nh + Nd 



Ljr, 



(47) 



If we defined the number of empty lattice sites as Nh — 
L — N + Nd we see that ^'^^■^ /L can be written as 



Nd 



L 



LNh + Nd L 



(44) 



In addition to the interesting behavior at half-filling, we 
also see from the above equation that if the system is not 
at half-filling but doubly occupied sites are allowed then 
the flux felt by the c, excitations is not the usual (p/L 
but the effective flux ^""^-f /L. 

The corrections oi 0{t^ /U) to the energy (^ 
are given by 



and 



d'^Em 



= -z'E2. 



(48) 



Away from half-filling, z is finite and /U corrections 
will produce no qualita tive changes in the results we have 
obtained in subsection II B . We have checked that the ef- 
fect oi t^/U corrections is to increase the value of D{T) 
over that of the U = oo case. This behavior is in agree- 
ment with the physical picture that transport of charge 
is easier when electronic correlations are reduced. 



oj-2 



(45) 



In the limit ?7 3> i it is simple to compute the ground- 
state energy as a function of for a system of any size 
L, which reads 



E = 



, , ,,,sin(7rAr/i) St'^NB'o 



IQt^B, 



ULN 



rsin'(0/i 



sin(7r/L) UL 
sm^{TTN/L) 



8fBf^ sm{2TTN/L) 
Ul sin(27r/L) 



sin2(7r/i) 
cos(20/L) , 



(46) 



with _Bq = A^ln(2)/2, in the thermodynamic limit. From 
equation ( ^ ) it is simple to obtain D{Q, L), for any band 
filling. 

It is instructive to consider Eq. ( |46[ ) in two different 
cases: (a) one electron off half-filling; (b) at half-filling. 
The first case gives for the charge stiffness D{0, L) = tS, 
where 6 — {L ~ N)/L = 1/i is the dopping; for the 
second case we obtain D{Q, L) — 0. Both these results 
agree with those of Stafford, Millis and Shastryll3. In the 




-0.1 0.1 
<I>/2jt 



FIG. 4. The B - iVc,iln2 term as function of the flux 
(j)/{2Ti) (0 is in units of the flux quantum 2%), for Nc,-i ~ 10 
(solid line) and Nc,i = 40 (dashed line). It is clear that 
dB/d<j) is zero at = 0, and, therefore, the corresponding 
current matrix element jm is zero. 
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At half-filling z = and, therefore, the current and 
the curvature of levels are, up to order t, zero for all the 
eigenstates of the model, leading to a zero value of D{TX. 
This result agrees with that of Fujimoto and KawakamiO 
if we take U — oo in their result. 

If, for electronic densities equal to one, t'^/U correc- 
tions are included, the flux dependence of the energy 
levels is due to the rapidities R'^jj only. The flux de- 
pendence of i?^^ j controls the behavior of with (p/L. 
This factor is the only responsible for the flux dependence 
of the energy eigenstates, as can be seen from Eq. (p5|). 

To simplify our study, we consider only the c, 1 string 
excitations. From Eqs. (|lj) and (jl^) it is clear that 
the creation of Nc^i charge-string excitations such that 
= Ni =Const introduces the constraints 

AN, = N,- TVf = -2A7V,,i = -27Ve,i 

TVf^iV, iVfo-^i- (49) 

where ANa^j represents the variation of the numbers 
Na^.y relatively to the ground state (only spinon exci- 
tations have been considered) and -/Vf *g is the number of 
the occupied quantum numbers /"'° in the ground state. 

From Eq.p5|) we see that the spin degrees of freedom 
do not contribute to the flux dependence of the energy 
levels since they do not depend in any way on 4>/L, up to 
order t'^ /U. All the flux dependence of the energy levels 
is in the term 

with B depending on the eigenstate m considered. 

As we have stressed, up to order P/U the charge 
gapped excitations are mathematically equivalent (apart 
from the gap relative to the c, 0-excitations) to an 
isotropic Heisenberg model in a lattice of size (L — Nc). 
Combining this similarity with Eq. (^9|), the sum B 
reads, for large iV^^i, B = iVt,^iln2, at (/) = 0. At finite 
4>, the sum B depends on </), but presents a maximum for 
(f> — 0, as can be seen from Figure ^. If we take, as an 
example, the case L — TVc = 2 it simple to see that Eqs. 
( p^ ) and (^) do reproduce the flux dependence of the 
energy of a two-site Heisenberg chain. 



III. ALGEBRAIC SOLUTION OF THE 
HUBBARD CHAIN WITH TWISTED 
BOUNDARY CONDITIONS 

While the knowledge of the low-energy eigenstates is 
enough if we are interested in the low-energy or temper- 
ature behavior of the Hubbard model, when calculating 
the conductivity of the Hubbard model at half filling, 
states with energies of the order of U become relevant. 
In this section, we obtain all eigenstates and eigenvalues 



of the Hubbard model with a magnetic flux, when t/U is 
a small parameter ptakpg an alternative approach to that 
of the BA methocB^ll§E3, used in the previous section. 

The eigenvalues of the Hubbard model, when t = 0, 
are given by En^ = N^-U and all states with a given 
number Nd of doubly occupied sites are degenerate. Us- 
ing the identity c.^r = c,cr[{l - rii^) + n^af^, where c.,,o- 
is the fcrmion annihilation operator on site i with spin 
a and ni ,j — cl^Ci^^, the Hubbard model, Eq. (^, with 
t 7^ can be rewritten a E 

H = fh + fd+ {fhd + fdh) + UY, n,,Tna (51) 
where 

2\ = -t ^ (1 - nig)c\^Cj„{l - njff), (52) 

Td = t ^ ^ nigC^^CjaTljg , (53) 

Tdh = X! "»*'^Lcjt(1 - njs), (54) 

with Tdh — Tj^d ^^'i ^ ~ Here, Td and Th describe 
the movement of double occupations and holes, respec- 
tively, but these terms and the Hubbard onsite interac- 
tion do not change the number of double occupations of 
a state. An applied magnetic flux is easily introduced as 
a twisted boundary condition. In the limit of U/t oo, 
the first correction to the t = spectrum is obtained by 
diagonalizing the model given by Eq. ( pl| ) without the 
mixing term Thd + Tdh within each subspace Nd- Such 
diagonalization is described below. 



A. Exact diagonalization using slave-fermions: U ^ t 
results 

Let us consider a ring of L sites (labeled clockwise) 
with L even and a general state with Nh holes and Nd 
double occupancies {Ns singly occupied sites) with posi- 
tions defined by the ordered sets {h} and {d} and with 
the spin configuration of the singly occupied sites given 
by the set {a}: 



{h} 






{d} 


= {di,.. 






= {cti, . . 


.,aL-Nh-Nj 









|{M;M;W) = n4T4ani-J0) (55) 
1=1 j=i 
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with bj = j + nj{{h}, {d}) where nj{{h}, {d}) is the num- 
ber of holes plus double occupancies to the left of site 
j. The Hamiltonian applied to this state exchanges a 
hole or a double occupancy with a spin, without chang- 
ing the spin configuration {a} except for hoppings at the 
boundaries where the final spin configuration is a cyclic 
permutation of the original one. ___ 
We now introduce the slave-fermion representationEHo'Cj 
for the fcrmion operators, Ci^^ — elSia + crSj _^di, where 
ej and dj are operators satisfying fermionic commuta- 
tion relations ({e|, d|} = 0) and Sia is a spinon operator 
satisfying bosonic commutation relations with the re- 
striction ejci + d\di + Sj^Sia — 1. Let us define the 
ordered set {a} — {h} U {d} = {ai, . . . , ajv^-i-Ard} and 
a set of fictitious spins {i^} = {i^i, . . . , j^Ar^+TVd} which 
indicates whether we have a hole or a double occupancy 
at Oi, that is, 

= el^; C^.^ = dl^. 

These operators obey fermionic commutation relations. 
Our general state in the slave-fermion representation, us- 
ing the previous definition, becomes 

|{a};M;M> = n<., 11 ^I.-JO)./ 

3=1 k=l 

X (-l)Er=i(^'-i' (56) 

H^Nd-U^Y. t^-S\^l,aS^..Cl^C,+ l,u + h.C. (57) 



and the Hamiltonian is 

L 



with tiy = = V - t for i ^ L and Ilv = e^'^''^ti^. A i/-spin 
dependence of the hopping integral just reflects that the 
hole band is shifted by tt relatively to the electron band. 
Note that with this representation, the boundary condi- 
tions also become i^-spin dependent, since a hopping of a 
hole from site L to 1 implies an opposite hopping of an 
electron from 1 to L. As we shall see, this will lead to an 
anomalous flux dependence of the U ^ t eigenstates. 
Let us define 

\a'i, . . . , ajVfc+Afd' '^'n ■ ■ ■ 1 ^L-Nu-Na) ~ ^^^) 

In order to simplify this Hamiltonian, we construct Bloch 
states in the cyclic permutations of the spin configuration 



\{h};{d};as,qs 
1 



(59) 



'|{/i};{d};ai_ 



where = n(27r/rQ,^) with n = 0, • • • , r^^ — 1 and r^^ is 
the period of the spin configuration. The Us index labels 



the different (not obtainable from any other by cyclic 
permutations) spin configurations with period Tq,^ . Then 



(60) 



(-1) \{a'};{iy'};a,,qs), i^L 
(-l)e»?= |{a'};{z/'};a„g,), i^L 



For the state given in Eq. (|56[), the successive applica- 
tion of Th and Td does not change the sequence of these 
fictitious spins {z^}, but permutates them cyclicly as in 
the case of the real spins {a}. Given a {as, (Zs} subspace, 
the Hamiltonian can be written as: 



+h.c.] +Nd-U 



7^ L, tLv 



(61) 



and 



with til, — 

Let us consider first the zero fiux case, </> = 0. In that 
case, the model above is the one-band {U — oo) Hub- 
bard model with twisted boundary conditions, except for 
the spin dependence of the hopping integral. In order to 
make tii = we make the following gauge transforma- 
tion: 



Then ti^ t, i ^ L and = e^'^^'t. Now, the gauge 
transformation 



leads to tji, te^^ , Vj, z^. The solution of this model 
is known (this is the U — oo Hubbard chain enclosing a 
magnetic flux (l)oqs/n), but now our up and down spins 
indicate whether the charge carrier is a hole or a dou- 
ble occupancy. The spin-charge factorization which we 
know characterizes the eigenf unctions of the LL — oo 
Hubbard model (we refer the readers to Ref.EZl for a 
detailed solution of this model) translates to a further 
decoupling of the part of the wavefunction describing 
the charge degrees of freedom and consequently, the 
eigenstates are factorized in the form: | charge carriers 
positions) (8) \ hole/ double occupancy configuration) ® \ spin 
configuration) (see Fig. ^). That is, the eigenstates of 
our model can be mapped onto a system of fermions 
on a chain with L sites and two squeezed chains, one 
of holes and double occupancies and the other of spins. 
The eigenstates will be of the form 

\{k}) (g) \ac,qc) ® \cts,qs) 

where qc is the momentum of the Bloch state in the cyclic 
permutations of the {v} conflguration, qc = ni^-K jr^^ 
with n = 0, • • • , ra^ — 1 and r^^ is the period of the 
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holes/double occupancies configuration. The charge car- 
riers will be free fermions in a chain threaded by a ficti- 
tious magnetic flux cf) — [q^ — qc)4>o/TT generated by both 
the spin and hole/double occupancy configurations and 
the eigenvalues of H{qs) — Nd ■ U are given by 



E{ki,---, kNu+Na ) ^ cos 



L 



(62) 



with k, = (27r/L)7ij, = 0, • • • , L - M The fact that 
the hole/double occupancy configuration also generates 
a fictitious flux implies that the lowest eigenvalue of a 
Nd 7^ subspace may have an even weaker dependence 
on an applied flux than that of the ground state energyE3. 





® 9 



(8) 



FIG. 5. The high energy eigenstates of our model are 
mapped onto a system of free fermions on a chain threaded 
by a fictitious fiux cf> — {qs — qc)4'o / and two squeezed chains, 
one of holes and double occupancies and the other of spins. 



In order to solve the (/> 7^ case, we first repeat the 
gauge transformations Cj| Cj^e"*'^^ and Cj^, 
t ^t^j so that tjy te'^^ , Vj ^ L and tn, 



te* "£ e 



Now, we define a operator Q such that 



(63) 



that is, it does a circular permutation of the spin configu- 
ration. For the </> 7^ case, given a configuration of ficti- 
tious spins {ly}, it is obvious the eigenstates will be found 
in the subspace spanned by Q"'{v}, n — 0, • • • , r^^ — 1. 
The hoppings across the boundary do a cyclic permuta- 
tion with a factor e*" "^. We wish to construct now the 
state such that 



J2 «41M> =e^*' J2 «'Q1M> , (64) 



where — Qe"^ if =t and 



if I'l =i. 



b' will be the effective flux felt by the noninteracting 



fermions. That is, we want to construct the equivalent 
to the Bloch states in the cyclic permutations of {ct} for 
the (f) — case. 

This problem is equivalent to solving a one-particle 
tight-binding model for a chain of Tq,^ sites with hop- 
ping constant tj — te^^^''^ with the correspondence \i) = 
Q^~^\{i^}) . The total flux through this tight-binding 
chain is 



Nh-Nd 
'Nn + Nd' 



(65) 



The solution is obtained after a gauge transformation to 
make tj e*'^i/'"°=t. The gauge transformation depends 
on the i/-spin configuration, but the tight binding eigen- 
values only depend on the total flux. The eigenstates 
will be Bloch states \ac,qc) (in the cyclic permutations) 
in the way defined for (j) — with qc = n{2TT/ra^) with 
?i = 0, • • • , — 1 and eigenvalues given by 



e{qc) = 2i cos ( qc 



This implies 



\ac,qc 



(66) 



(67) 



and therefore, repeating the (j) — steps, the Hamiltonian 
can diagonalized and the eigenvalues of H{qs) — Nd • U 
are given by 



E{ki,- ■ ■ ,kN^+Nj ^ 2t ^ cos 



L 

1 Nh-Nd ^ 
' LNd + Nh 



(68) 



with ki ~ {2Tr/L)ni, ni = 0, — 1. Summing up, 

we have derived an effective Hamiltonian for our spinless 
and chargeless particles 



i/(q„g,)=te-'*^4A+i+h.c. 

i 

with 

and and ^'^^ ^ given by 

^q,-qc, 

and 



Nh~Nd 
Nd + Nh^ 



(69) 

(70) 
(71) 

(72) 



The Hamiltonian describes free fermions in a chain 
threaded by a fictitious magnetic flux 4>o{gs — <1c)/t^ 
(see Fig. H). The lowest energy state is obtained for 
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qs — Qc = and Qs — Qc = TT for an odd and even num- 
ber of holes plus double occupancies respectively in the 
case of zero external flux. The flux given by Eq. ( |70| ) 
has exactly the same form as obtained previously. The 
^stat ^Qj-jj^ corresponds to the sums over the Ij''' and Ij'^ 

numbers in Eq. (p5|), and the term ^'^f-f is the same as 
in Eq. (Q). It quite clear from this approach that the 
presence of doubly occupied sites is responsible for the 
renormalization of the flux felt by the holons. 



B. Transport of charge 

In the limit t ^ U, an effective Hamiltonian to sec- 
ond order in t can be obtained for each Nd suibspace by 
a canonical transformation Heff = e^^ He~^^c3 leading 

to the additional term {fhdfdh - fdhfhd)/U, {Nd ^ 0). 
This term will give a correction of order t'^ /U and lift the 
degeneracy of states \{k};ac,qc]C(s,cis) in what concerns 
the spin and charge configurations \ac, Qc] cxs, Is) as in the 
case of Nd — subspace where a similar term-leads to a 
Heisenberg interaction plus a three-site termO'EZI. 

The Hubbard gap, that is, the energy difference be- 
tween the highest eigenvalue with Nd double occupancies 
and the lowest eigenvalue with Nd + 1 double occupancies 
for a system with N electrons, is easily obtained from Eq. 



(73) 



with (5 = for Nh + Nd odd and 6 = tt/L for Nh + Nd 
even. For a half-filled system, N^ — Nd and the gap 
grows with Nd- Eq. is the same as Eq. (|l]), if t^/JJ 
corrections are neglected. 

The Hubbard Hamiltonian can be written as a sum 
of a diagonal term in the number of double occupancies 
and a mixingr-term which has little effect in the strong 
coupling limitcZI. Similarly, for the current operator. 



■C!+l,cr 



(74) 



1 

71 



('^i+l,T'5'i,x 



(77) 



(78) 



Note that in Jo, no i^-spin dependence appears in the 
hopping integral because the term in this current oper- 
ator that hops a double occupancy from i to i -I- 1 has 
a minus sign while the term that hops a hole from i to 
i -V \ does not. The Jo current operator is now simply 
the double occupancy current minus the hole current. 

One should note that when determining the optical 
conductivity, one has to calculate matrix elements of the 
current operator between states with Nd and Nd + 1 dou- 
ble occupancies in order to obtain the upper Hubbard 
band part of the optical conductivity. Obviously, we can 
replace J by J\ in this case. The low frequency region is 
given by matrix elements of the current operator between 
states with the same number of double occupancies and 
in this case, we replace J by Jo. 

The expectation value of the current operator for our 
zero flux eigenstates can be easily obtained if we rewrite 
the current operator as 



Jo 



(2Cl) ■ + h.c. (79) 



where Cf = {C\^Ci^-\ — C\^Ci^\)l2 and with ti^ given 
as in Eq. This term just means that a hole current 
implies an electron current in the opposite direction, but 
a double occupancy current is an electron current. The 
second term is of the same form as that in Eq. ( p7|) and 
applied to an eigenstate, it would multiply it by a phase 
if the CI term was not present. One should note that if 
we consider the squeezed i/-spin chain, will 
give the value of the spin in site m, but (ckc, (?c|Cf„|ac, 9c) 
is independent of m since the state |q;c,<7c) is invariant 
under a t^-spin translation. We can therefore replace it 
by a sum over all sites in the reduced ^'-spin chain, di- 
vided by the total number of sites of this squeezed chain, 
Nd + Nh ■ And obviously. 



we can write J ~ Jo + Ji where Jo commutes with 
H{(j) — 0) and Ji is the mixing term. These two op- 
erators in the slave-fermion representation are given by 

Jo = it {S\+l.(ySi,(yC\.vCi+l,i' ~ Sl^^Si+i^aCl^i ^^Ci^u) 



and 



(75) 



Ji = i2tY,{Cl+i.A+i, - C.+i,.5|+i,J (76) 



where 



1 



Nd + Nh 



Nh~Nd 
Nh + Nd 



(80) 



Clearly, the average value of Ji is zero and so, using the 
above expression for Jo, we obtain 



({fc}; Uc, qc, as,qs\J\{k}; a^, qc, a^, qs) = 



(81) 



Nh + Nd ^ 



sm 



The equation derived above for the mean value of the cur- 
rent operator is the same as Eq. (^7|), derived from the 
BA, it leads to D{T) = 0, at half fiUing, and to D{T) ^ 
otherwise. 
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This result and the ones described previously suggest 
a simple picture for the U ^ t Hubbard model, that of a 
tight-binding model of + fermions with an average 
charge {Nh — Nd)/{Nh + Nd) moving in a averaged spin 
background (5^) = 1/2 • {N^ - Ni)/{[N^ + N^). 



IV. COMPARISON WITH AN IDEAL SEMI 
CONDUCTOR 

The simplest one-dimensional model-Hamiltonian for 
an ideal semi-conductor is the dimerized tight-binding 
model, which, after diagonalization, reads (the spin in- 
dex is omitted)Ea 



k,l3 



with 



2tit2 cos(2fca) 



(82) 



(83) 



and ^1,2 arc the hopping integrals, a the lattice spacing, 
f3 the band index, and k the momentum, belonging to the 
Brillouin zone [— 7r/(2a), 7r/(2a)]. The zero-temperature 
charge stiffness is easily computed, and is given by 



D{0) = 



4Ut 



11^2 



sm(kpa) cos{kpa) 



(84) 



with kp the Fermi momentum. At zero temperature, 
D{0) is zero at half filling, the contributions to D{0) be- 
ing only from the lower band (3 = —I. At finite tem- 
perature there are contributions to D{T) from particles 
in both bands. However, and at odds with the Hubbard 
model, finite occupancies of the upper band (/3 = 1) does 
not renormalizc the flux felt by the holes that remain in 
the lower band. For the dimerized model we obtain a 
finite charge stiffness at half filling, given by 



2 /"f/(2<i) 
D{T)^- dkf{k,T)d{k), 

J-Tr/(2a) 



with d{k) given by 

d{k) = - 



£;4(/c) 



2\2 



E^{k) 



(85) 



(86) 



and f{k,T) = (l-Hexp(i;(fc)/r))-i (the chemical poten- 
tial is zero). We plot D{T) in Figure |^. From this Figure 
we see that the maximum of D{T) is controlled by the 
energy gap A — 2\ti—t2\. In the Hubbard model had not 
the effective flux felt by the holons been renormalized to 
zero at half filling, and we would have obtained for that 
model a physical behavior similar to that in Figure ^. At 
T — and n = 1, we can still view the Hubbard model 
as a hand insulator, since the Mott-Hubbard transition 
has its origin in the gap between the bands of c, and of 
the c, 7 excitations. 




FIG. 6. Charge stiffness for the dimerized model, at 
half-filling. The parameters are a = 1, t\ = 1, t2 = 2. The 
temperature axis is scaled by the gap A. 



V. CONCLUSIONS 



In this paper we have studied the behavior of the finite- 
temperature charge stiffness for Hubbard rings of finite- 
size L. 

For U — oo and n < 1, we have seen that D{T,L) is 
finite, decreasing as the temperature increases. In finite 
rings there is a difference between free spinless fermions 
and the c, excitations of the Hubbard model, due to the 
statistical flux (/)^*'**. This statistical flux is also respon- 
sible for the sensitivity of D2{T, L) - the contribution to 
D{T,L) from the square of the matrix elements of the 
current operator ~ to the inclusion of all the spin exci- 
tations. At high temperatures we have seen that this 
sensitivity vanishes. We have also verified that D{T,L), 
as calculated from the curvature of levels Dm, is not very 
sensitive to the inclusion of only some of the spin excita- 
tions, even at low temperatures. At half-filling, the only 
available state, apart from the massive spin degeneracy, 
is the ground state which leads to D{T) — 0. Our Jesuits 
agree with the statement of Zotos and Prelovsektj that 
D{T) should be finite if ^(O) is also finite. 

For J7 3> t, we have derived asymptotic equations for 
the mean value of the current operator and for the cur- 
vature of levels valid for any band filling. This derivation 
is specially important at half-fiUing, where some conjec- 
tures by Zotos and Prelovsektj indicate that D{T, oo) 
should be zero. These authors related this odd behav- 
ior with the integrability of the model. We have shown, 
using two independent methods, that if doubly occupied 
sites are allowed the flux felt by the low-energy charge 
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excitations - holons - is renormalized to zero, at n = 1. 
For this density, only the gapped charge-excitations feel 
the flux piercing the ring, but the current transported 
by these excitations is zero, and therefore, D{T) — 0. 
We stress that his result does not coiitradict what was 
obtained by Fujimoto and KawakamiO, since they also 
obtain D{T) = in the limit we are considering. Nev- 
ertheless, our analysis shows that if D{T) = is to be 
finite, for n = I, the prefactor of the exponential of the 
Mott-Hubbard gap must be, at least, of the order of t'^/U, 
and not oij-the order one, as it is found by Fujimoto and 
KawakamS. The physical interpretation of this effect 
is quite transparent in our approach: the higher energy 
excitations are renormalizing to zero (for n = 1) the flux 
felt by the low-energy excitations. In addition it is also 
clear from our approach that the highest occupied band 
is always completely filled and has a maximum at (/) — 0. 
This effect, together with the renormalization of the flux 
leads to a zero value of the current, in any eigenstate 
of the model. Another way of looking at our results is 
that the U :S> t Hubbard model, can be thought as a 
tight-binding model of + Nfi fermions with an aver- 
age charge {Nh — Nd)/{Nii + Nd) moving in a averaged 
spin background (S*^) = 1/2 • {N^ - Ni)/{{N^ + Ni). 

If the system is not at half filling and doubly occupied 
sites are allowed, the flux felt by the holons is renormal- 
ized away from the usual 0/L value by the presence of 
doubly occupied sites. 
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quantum numbers 


type of excitations 


gap 


rapidity 




holons 


gapless 


kj 


js,0 


spinous 


gaplcss 






charge strings 


gapr)ccl 




c 


hpiu stings 


gapped (7/ / 0) 





TABLE I. Types of excitations associated with the quantum rmnibcrs with a = c,s. The spin-string excitations are 

gapped only for finite magnetic field H. Both the holon and spinon excitations can be described by bosonization methods. The 
smaller gap associated with the charge string exciations is the Mott-Hubbard gap Amh- 



state Nc Nsfi Nc,j>o Ns,j>o 



G. S. 


N 


Ni 








Ex.0 


N 












N - 2'^No,-y 


Ni - 7^c,7 









N 


iVi-(7 + l)iV,,^ 








TABLE II. The numbers, Nc,-/, of the different kinds of excitations, in different classes of cigenstates. The notation is as 
follows: G. S., Ex.Q, Ex.c,-y, and Ex.s,-^ stand for the ground state, the low-energy eigenstates (no-strings), the eigenstates with 
Nc,j charge strings of length 7, and the eigenstates with Ns,-y spin string of length 7, respectively. 



states (J|)'"'" Nc (J^^.O)'nm jv^^^ (j-a,l^min jy^^^ ^js,2^min jy^^^ ^ 



Ex.0 


-11/2 


6 


^ .7 ' 
-1 


3 


^ .7 ' 




^ .7 ' 







Ex.s,i 


-5(-6) 


6 


-1 


1 





1 






— — — 

36 ' 36 


Ex.s,2 


-11/2 


6 













1 






TABLE III. The notation is as follows: Ex.o, Ex.s,i, and Ex.s,2 stand for the states with no spin strings, the eigenstates 

with one 7 = 1 spin-string, and the eigenstates with one 7 = 2 spin-string, respectively. The system size is L = 12, the number 
of electrons is Nf = Ni = 3, and U = 00. The number of occupied quantum numbers I^''' is given by the Na,^. The numbers 

of available /j^'^ is given by the set {(/"•T)'"*", (7«,7)min + 1^ . . . , (7«.7-)max| ^^^^ (ja.T^min ^ _(ja,7)max^ gxpgp^. foj. 

umklapp processes can break this symmetry. The column 01*"* indicates the possible statistical fluxes. 



states 




Nc 








iV.,0 




iV.,1 


^j-s,2-jrnin 


Ex.0 


-17/2 


6 






-1 


3 








Ex.s,i 


-8(-9) 


6 






-1 


1 





1 




Ex.s,2 


-17/2 


6 

















EX.(c,l;s,l) 


-8(-9) 


4 


-6 


1 









1 




Ex.c,i 


-17/2 


4 


-6 


1 


-1/2 


2 








Ex.c,i 


-17/2 


2 


-13/2 


2 





1 








Ex.c.i 






-7 


3 












Ex.c,2 


-8(-9) 


2 


-6 


1 





1 








Ex.c,3 






-6 


1 













TABLE IV. The notation is as follows: Ex.q, Ex.^^i, Ex.a,2, and Ex.(c,i;s,i) stand for the states with no spin strings, the 
eigenstates with one 7=1 spin-string, the eigenstates with one 7 = 2 spin-string, and the eigenstates with one 7=1 charge- 

and one 7 = 1 spin-strings, respectively. Ex.c,-, stands for states with only charge- strings. The system size is L = 18 and the 
number of electrons is Ni = Ni = 3. The number of occupied quantum numbers 7"'^ is given by the A^a,^. The numbers of 

available 1°'''^ is given by the set {(/"'T)-"*", (7«.7)mm + 1, . . . ^ (ja,7)max| g^^^ (ja,7)min ^ _(ja,7)max^ g^pgp^ foj. 

unklapp processes can break this symmetry. 
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